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^ ! Abstract 

In this work we develop a new propagation model for the Galactic cosmic rays based on the 
GALPROP code, including contributions from dark matter annihilation. Its predictions of the 
Galactic diffuse 7 ray spectra are compatible with the EGRET data in all sky regions. It also gives 
consistent results about the diffuse 7 ray longitude and latitude distributions. The B/C, 10 Be/ 9 Be, 
proton, electron and antiproton spectra are in agreement with cosmic ray measurements as well. In 
this model we have taken a universal proton spectrum throughout the Galaxy without introducing 
large fluctuation, considering the proton energy loss is negligible. The dark matter annihilation 
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signals are "boosted" after taking the contributions from subhalos into account. Another interesting 
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feature of this model is that it gives better description of the diffuse 7 rays when taking the source 
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distribution compatible with supernova remnants data, which is different from previous studies. 
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I. INTRODUCTION 



Cosmic ray (CR) propagation is a complex process involving diffusion by magnetic field, 
energy losses and spallation by interactions with the interstellar medium (ISM). Diffuse 
Galactic 7-rays are produced via the decay of neutral pion and kaon, which are generated 
by high energy cosmic nuclei interacting with interstellar gas, and via energetic electron 
inverse Compton (IC) scattering and bremsstrahlung. The 7 rays are not deflected by the 
magnetic field and the ISM is transparent to 7-rays below a few TeV pj. Therefore, the 
observation of the diffuse 7-ray spectra and distribution is a valuable diagnosis of the self- 
consistency of propagation models, the distribution of CR sources and the ISM. The Galactic 
diffuse 7 rays has been measured by EGRET 

as and exhibits an excess above ~ 1 GeV 
compared to prediction [3]. The theoretical calculations are based on a conventional CR 
model, whose nuclei and electron spectra in the whole Galaxy are taken to be the same as 
those observed locally. 



since it was first raised. 



The discrepancy has attracted much attention 
It may either indicate a non- ubiquitous proton or electron spectrum, or the existence of new 
exotic sources of diffuse 7-ray emission. 

Many efforts have been made to solve the "GeV excess" problem within the frame of 
CR physics, such as adopting different CR spectra , or assuming more important 

contribution to diffuse 7-rays from CR sources [7]. A brief review of these efforts is given in 
{f|. In that paper an "optimized" propagation model has been built by directly fitting the 
observed diffuse 7-ray spectrum. This "optimized" model introduces interstellar electron 
and proton intensities that are different from the local ones and reproduces all the CR 
observational data at the same time. Up to now, it seems to be the best model to explain 
the EGRET diffuse 7-ray data based on CR physics. However, this "optimized" model is 
fine tuned by adjusting the electron and proton injection spectra, while keeping the injection 
spectra of heavier nuclei unchanged, as in the conventional model, so that the B/C ratio is 
not upset. Furthermore a large scale proton spectrum different from the locally measured 
one might not be reasonable, since the proton diffusion time scale is much smaller than 
its energy loss time scale, which tends to result in a large scale universal proton spectrum 
within the Galaxy apart from some specific sources. Unlike protons, the electron spectrum 
may have large spatial fluctuation due to their fast energy losses from IC, bremsstrahlung, 
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ionization and the stochastic sources 111 ]. 



Another interesting solution, given by de Boer et al. Q, [hJ, is that the "GeV excess" 
is attributed to dark matter (DM) annihilation from the Galactic halo, where the DM 
candidate is the neutralino from the supersymmetry (SUSY). By fitting both the background 
spectrum shape from cosmic nuclei collisions and the signal spectrum shape from dark 
matter annihilation (DMA) they found the EGRET data could be well explained 8j. This 
suggestion is very interesting and impressive, due to the fact that in 180 independent sky 
regions 1 , all the discrepancies between data and the standard theoretical prediction can be 
well explained by a single spectrum from DMA with m x = 50 ~ 70 GeV. Furthermore, 
by fitting the spatial distribution of the diffuse 7-ray emission they reconstructed the DM 
profile, with two rings supplemented to the smooth halo. The ringstructure seems also 
necessary to explain the damping in the Milky Way rotation curve 8| and the gas flaring 
jlo| . However, the DMA solution to the "GeV excess" also meets a great challenge because 
of its prediction of the antiproton flux. In de Boer's model, this flux is more than one order 
of magnitude greater than data [3]. The overproduction of antiprotons comes from two 
factors: a universal "boost factor" ~ 100 of the diffuse 7-rays boosts the local antiproton 
flux by the same amount; the two rings introduced to account for the diffuse 7-ray flux 
enhance the antiproton flux greatly since they are near the solar system and are strong 
antiproton sources. In their work, de Boer et al. did not try to develop a propagation 
model. Instead they focused on reconstruction of the DM profile by fitting the EGRET 
data. They need a "boost factor" to enhance the contribution from DMA. The background 
contribution from pion decay is arbitrarily normalized in order to fit data best. 

In the present work we try to build a propagation model to explain the EGRET diffuse 
7-ray data based on both Strong's and de Boer's models while overcoming their difficulties. 
In our model the diffuse 7-ray comes from both CRs and DMA directly. On one hand we do 
not introduce a different interstellar proton spectrum from the local one; on the other our 
model gives consistent p flux even when including contribution from DMA. Furthermore we 
do not need the large "boost factor" to DMA or renormalization factor to CR contribution. 
Actually, the 7-ray flux from DMA is boosted by taking the subhalos into account. The 



1 The division at low energies may be problematic as the EGRET point spread function is about 6° and 
non-Gaussian at low energy. 
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: Galactic longitude 



131 ] . we have briefly 



diffuse 7-ray spectra at different sky regions and its profiles as functions o 
and latitude are well consistent with EGRET data. In a previous paper 
introduced our model. Full details are given in the present paper. 

The paper is organized as follows. We describe the calculation of the DMA contribution 
in Section II. In Section III, we focus on the conventional CR model. As underlined, it 
explains the EGRET data, but produces too large p flux. In Section IV, we present our new 
propagation model and its predictions for the diffuse 7-ray spectra and profiles, and the p 
flux. Finally we summarize and conclude in Section V. 



II. DARK MATTER ANNIHILATION 



Here, we calculate the diffuse 7-ray flux from DMA. The frame of minimal supersymmetric 
extension of the standard model (MSSM) is retained, where we assume that DM consists 
of the lightest neutralinos. A pair of neutralinos in the Galactic halo can annihilate into 
leptons, quarks and gauge bosons. Their decay products include a 7-ray continuum and 
thus contribute to the 7-rays diffuse emission produced by CR interaction with the ISM. 



A. Particle factor 



The flux of 7 rays from the neutralino annihilation is given by 

_ /cA (erf) dN /■ p 2 1 (av)dN 1 f 2 . . 2j JO _ 
W = U dE J/ dV A^ = ^kklE X * J P^ydrdV , (1) 

where (av) is the averaged neutralino annihilation cross section times relative velocity, ^ 
is the differential flux in a single annihilation, m x is the mass of neutralino, d is the dis- 
tance between the detector and the 7-ray source, and p(r) is the spherically-averaged DM 
distribution, determined by numerical simulation or by observations. The flux in Eq. ([T|) 
is determined by two independent factors, the first one only depending on the DM particle 
nature (mass, strength of interaction and so on), the second one depending on the DM distri- 
bution only. The first factor is denoted as "particle factor" and the second as "astrophysics 
factor" . 

The "particle factor" is calculated by doing a random scan in the SUSY parameter space 
and choosing models which could satisfy all the collider and cosmological constraints. How- 
ever, there are more than one hundred free SUSY breaking parameters, even for the R-parity 
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FIG. 1: g^**— >-y) as a function of neutralino mass. Thresholds for p and 7-ray kinetic energies are 
taken as 0.5 GeV. 



conservative MSSM. A general practice in phenomenological studies is to assume some sim- 
ple relations between the parameters. Following the assumptions in DarkSUSY [14] we take 
seven free parameters during the calculation of DM production and annihilation, i.e., the 
higgsino mass parameter /i, the wino mass parameter M 2 , the mass of the CP-odd Higgs 
boson ttia, the ratio of the Higgs Vacuum expectation values tan/?, the scalar fermion mass 
parameter mj, the trilinear soft breaking parameter A t and A b . All the sfermions are taken 
with a common soft-breaking mass parameter m^; all trilinear parameters are zero except 
those of the third family; the bino and wino have the mass relation, M\ = 5/3 tan 2 $w M 2 , 
coming from the unification of the gaugino mass at the grand unification scale. 

A random scan in the 7-dimensional parameter space of MSSM is performed using the 
package DarkSUSY [yj]. We assume these parameters to range as follows: 50 GeV < 
M 2 , M A , mj < 10 TeV, 1.1 < tan/? < 61, -3m§ < A t , A b < 3m 9 ~, sign(/i) = ±1. We 
find the 7-ray spectrum with m x = 40 — 50 GeV (after being added to the background 7-ray 
spectrum) can fit the EGRET data very well, which is consistent with the result of de Boer 

n 

et al. [8J. However, the branching ratio to 7 rays varies by a few order of magnitude for 
different models, and only models with large branching ratio into 7-rays can account for the 
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"GeV excess" . We also calculate the branching ratio into antiprotons since the p flux is a 
sensitive test of the DM model. 

In Fig. [T] we give the ratio of neutralino annihilation into p and 7-rays as a function of 
neutralino mass. We find g^**Z^j ~ ^ — jq above the threshold energy E t h = 0.5 GeV for 
m x = 40 — 500 GeV. For light neutralinos, the variation of the ratio in the parameter space 
is at most as large as a factor of 2. For heavy neutralinos the variation is very small, and 
the reason is that p and 7-rays are coming from the same final states, thus they are closely 
related. Below, we choose m x = 48.8 GeV to explain the EGRET data which predicts 
Qh 2 = 0.09 and Ijlfep^ ~ 0.055 for the threshold energy E th = 0.5 GeV. Actually for m x 
between 40 ~ 50 GeV the annihilated 7-ray spectra are almost identical. A few models also 
show similar branching ratio as in our chosen model (e.g., for m x = 41.13,41.21 GeV). 



B. Astrophysics factor 



2 

The astrophysics factor determining the annihilation fluxes is defined as $ astro = J v ^dV 
with d the distance to the source of 7-ray emission, p the density profile of DM, and V the 
volume in which annihilation taking place. In the de Boer model, the authors adopt a cored 
Galactic halo model with two DM rings, boosting the 7-ray flux by a factor of the order 



100. In our work we take a cuspy NFW 15] (or Moore 161 ]) profile. Especially we take 



the contribution from subhalos into account to enhance the annihilation signals. 

High resolution simulations of cosmological structure evolution reveal that in the cold 
DM scenario the structures form hierarchically and a large number of substructures survive 
in the galactic halos. A fraction of about 10% of the total halo mass may have survived 
tidal disruption and appear as distinct and self-bound substructures inside the virialized 
host halos. The existence of substructures will enhance the annihilation rate greatly by 
enhancing the astrophysics factor in Eq. (CQ). 

The mass function and spatial function of subhalos are given by N-body simulations. A 
simple analytical fit to simulations allows to write the probability of a subhalo with mass m 



appearing at position r as 




-1 



n(m } r)=n ^—) ( I ■■■( — ) j ■ (2) 
where n is the normalization factor so that about 10% of the halo mass is enclosed in 



subhalos, M V i r « 1.0 x 10 12 M Q is the Galactic halo mass, Th = 0.14r v i r m 29 kpc is the 
core radius for the distribution of subhalos. The result given above agrees well with that of 
another recent simulation by Gao et al. 
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The DM density profile within each subhalo is taken as the NFW [15j, Moore [16] and 
a cuspier form |19j as p = j—r ym^T^^y with 7 = 1.7. The last form is favored by the 
simulation conducted by Reed et al. [20|, which gives that 7 = 1.4 — 0.08 log(M/M*), for 
the subhalo mass of 0.01M* ~ 1000M* with a large scattering, increasing as the subha 
mass decreases. Small halos with large 7 = 1.5 — 2 are also found by Diemand et al. 
We take 7 = 1.7 for the whole range of subhalo masses as a simple approximation. 

We calculate the concentration parameter c v by adopting the semi-analytic model of 
Bullock et al. [22], which describes it as a function of the viral mass and redshift. We 
adopt the mean c v — m su b relation at redshift zero. The scale radius is then determined 
by r™' w = r v /c v , r™ oore = r"^ w /0.63 and r] = r™* w /(2 — 7) for the three density profiles 
respectively. Another factor determining the 7-ray flux is the core radius, r core , within which 
the DM density should be kept constant due to the balance between the annihilation rate and 
the rate of DM particles falling into this region 23[ . The core radius r core is approximately 
10~ 8 ~ 10 -7 kpc for 7 = 1.7 and about 10" 9 ~ 10~ 8 kpc for the Moore profile. 

Along a direction (6,(f)), the subhalos contribute to the "astrophysical factor" $ astro = 
J t o s ^ sub dN su b(6, </>), where $ sub = J v ^47rr 2 dr is the astrophysical factor for a single 
subhalo, and N su b is the number density of subhalos. When we calculate the integration 
along the line-of-sight starting from the Sun, we get the Jacobian determinant as — , with r' 
the distance from the subhalo to the Sun. The minimal subhalos can be as light as 1O _6 M 
as shown by the recent simulation conducted by Diemand et al. 2l|, while the maximal 
mass of substructures is taken to be Q.01M V [24|. The tidal effects are taken into account 
based on the "tidal approximation" [24|, where subhalos are disrupted near the Galactic 
center (GC). The total signal flux comes from the annihilation in the subhalos and in the 
smooth component. 

In Fig. [21 we show the factor $ astro from the smooth component, the subhalos and the 
total contribution as a function of the direction to the GC. The Q astTO from subhalos is 
almost isotropic in all directions as the Sun is near the GC compared with the viral radius 
r V i r of the Galactic halo. The largest enhancement for 7 = 1.7 subhalos is observed at large 
angles and can reach 2 orders of magnitude. This enhancement depends on the value of 
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FIG. 2: The astrophysical factor Q astro (in unit of (GeV / cm 3 ) 2 kpc Sr 1 ) from different directions. 
The almost horizontal lines correspond to the contributions from subhalos only. 

?"core, while for the Moore profile the enhancement is about one order of magnitude, and for 
for the NFW profile only about 3 times larger. The $ astro for Moore and NFW profiles is 
not sensitive to r core 24(. We also notice that near the GC there is no enhancement. This 
is actually a very important difference from the model given by de Boer where the "boost 
factor" is universal [8|. Given the factor $ astro and the SUSY model we can predict the 
7-ray flux from neutralino annihilation. 

In the next section we give the background diffuse 7-rays from CR interactions with the 
ISM. 

III. THE CONVENTIONAL MODEL 



The optimized model of Strong et al. reproduces the diffuse 7 rays assuming interstellar 
proton and electron spectra different from those locally measured. However, the required 
fluctuation of the proton spectrum may be not realistic. The works of de Boer et al. have 
strongly indicated that DMA may account for the diffuse 7 ray excess [8]. Therefore our 
first attempt is to build a propagation model including contribution from DMA based on 
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the conventional CR model assuming universal CR spectra. 

We adopt the GALPROP package 0, Q to calculate the propagation of CRs and produc- 
tion of Galactic background diffuse 7 rays. It was shown that an explanation of the EGRET 
data in all sky directions in the conventional model is not an easy task, even including the 
contribution from DMA — considering that DMA provides only a single extra spectrum from 
neutralino annihilation. To give the best fit to the EGRET data, de Boer et al. [8] have to 
introduce arbitrary renormalization factors in different directions for the 7 ray spectra given 
in the conventional model 6J. Since we are trying here to build a propagation model, we will 
not introduce any renormalization factors, but rather explain the 7 ray data by adjusting 
the propagation parameters after adding the contribution from DMA. 

After a lot of tests, we came to the following propagation parameters. The scale height 
of the propagation halo Zh takes the same value 4kpc as that taken by Strong et al. in their 
conventional and optimized models 6|. The nuclei injection spectra share the same power 
law form in rigidity, and nuclei up to Z = 28 and all relevant isotopes are included. The 
CR injection spectra are given in Table [B We adopt the diffusive reacceleration propagation 
model. The spatial diffusion coefficient is given as a function of rigidity in the form 

D(p)=(3D (p/p ) 5 , (3) 

where (3 = v/c, D = 5.4 x 10 28 cm 2 s -1 , p = 4GV, and S = 0.34. The Alfven speed to 
describe the reacceleration process is Va = 34 kms -1 . These propagation parameters well 
describe the observed B/C ratio, the 10 Be/ 9 Be ratio, and the local measured proton and 
electron spectra, as shown in Fig. [31 

A major uncertainty in calculating the diffuse Galactic 7-ray emission is the distribution 
of molecular hydrogen, since the derivation of H 2 density from the CO data is problematic 
38j . The CO-to-H2 conversion factor Xco from COBE/DIRBE studies given by Sodroski 



et al. [39J] is about 2 — 5 times greater than the value given by Boselli et al. [40J in different 
Galactocentric radius. The later is based on the measurement of Galactic metallicity gradient 
and the inverse dependence of Xqo on metallicity [38J. The value of Xco is then normalized 



to the 7-ray data 38] . Strong et al. have derived the Xco by fitting the EGRET diffuse 7-ray 



data directly [4l|. A constant A co = (1.9 ± 0.2) x 10 20 cm" 2 /(K km s" 1 ) for E 7 = 0.1 - 10 
GeV was given 4l|. However, observations of particular local clouds yield lower values 
X co = 0.9 - 1.65 x 10 20 cm- 2 /(K km s" 1 ) [40]. Since the fit by Strong et al. to the EGRET 
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FIG. 3: Cosmic ray results in our conventional model. Lower and upper curves for B/C ratio 
correspond to solar modulated and local interstellar (LIS) values respectively. For 10 Be/ 9 Be ratio, 
the flat curve near 0.1 GeV/nucleon corresponds to the modulated flux and the other one is LIS. 



On the contrary, the lower curves for proton and electron spectra are the modu 
the upper are LIS. B/C data are from AC 
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TABLE I: Cosmic ray injection spectrum parameters. 
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FIG. 4: Cosmic ray source distribution along the Galactocentric radius R. The source distribution 
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is in arbitrary units. Dash-dotted line: based on pulsars 
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and in our conventional model[13j; solid line: adopted in our new model. Vertical bars are SNR 



data points 



43( | . These distributions are normalized at R=8.5kpc. 



data in 4jJ assumes only the background contributions, we expect that it gives a larger Xco 
than is the case with the new DMA component. We find that a smaller Xco ~ 1-0 x 10 20 
molecules cm~ 2 /(K km s" 1 ) can give a much better description of the EGRET data below 
1 GeV. After taking DMA into account, the global fitting is greatly improved. We assume 
for Xqq a constant value independent of the radius R. As shown in Ref. j38], the simple 
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form of Xco is compensated by an appropriate form of the CR sources. We have taken the 
radial distribution of CR sources in the form of {j^) a exp(— /3( - H ^"J ?Q - > ), with a = 0.5, (3 = 1.0, 
Rq = 8.5 kpc, and limiting the sources within R ma x = 15 kpc, which are adjusted to best 
describe the diffuse 7-ray spectrum. The source distribution is shown in Fig. H] (dashed line). 

Following Strong et al. 6j, we divide the whole sky into six regions: region A corresponds 
to the "inner radian", region B is the Galactic plane excluding the inner radian, region C 
is the "outer Galaxy" , regions D and E cover higher latitudes at all longitudes, and region 
F describes the "Galactic poles". The calculated diffuse 7 spectra in the six different sky 
regions are given in Fig. [5j The diffuse 7 ray background includes contributions from tt° 
decays produced by nuclei collisions, IC scattering off the interstellar radiation field (ISRF), 
bremsstrahlung by electrons, and the isotropic extragalactic 7 ray background (EGRB). For 
the regions A, B, C and D, i.e. the Galactic plane and the intermediate latitude regions, the 
it decay contribution is dominant. At high latitude regions E and F the EGRB becomes 
more important. To account for the GeV excess the peak of the DMA 7 ray spectrum 
has to have similar magnitude as the background. After adding the diffuse 7 ray emis- 
sion from DMA to the Galactic background, we obtain a perfect agreement with EGRET 
measurements for all sky directions. 

The result in Fig. [5] is really a success, considering that we simply add the two kinds of 
diffuse 7 ray contributions together, without introducing any arbitrary normalization for the 
background 7 rays or "boost factors" for the DMA contribution. It should be noted that 
including the enhancement by subhalos dose not exclude the ring-like structures proposed 
by de Boer [8|, |10J|. That is natural since taking the subhalos into account only enhances 
the signals coming from the smooth component, but does not mimic the ring-like structure, 



which can fit the EGRET data at different directions 



A. Flux of p 



a. 



We have seen in Fig. [T] that the branching ratios to 7 rays and to p are closely related. 
Therefore the p flux is a sensitive test of the DMA scenario to solve the "GeV excess" 
problem. In [12| , Bergstrom et al. have claimed that the de Boer model has been ruled out 
due to the overproduction of p flux. In their work they adopted a simple propagation model 
similar to that adopted in DarkSUSY 14j. We now check the p flux in our propagation 
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FIG. 5: Spectra of diffuse 7-rays for different sky regions (top row, regions A, B, middle C, D, 
bottom E, F). The model components are 7r° decay, inverse Compton, bremsstrahlung, EGRB and 
DMA (dark red curve). 
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model. 

We first calculate the source term of p produced by neutralino annihilation adopting the 
same SUSY model used for the diffuse 7 ray data, 



$p(r,E) 



2m^ 



-(p(rY) 



(4) 



where <f>(E) is the differential p flux at energy E due to a single annihilation, and 

(p(r) 2 ) = p\ moot h + (Psub) J 



(5) 



are contributions from the smooth DM component and subhalos. The contribution from 
subhalos is given by 



(p( r )Lb) = / n(m,r) ( / p 2 dV J dm, (6) 
where n(m, r) is the number density of subhalos with mass m at radius r, and J p 2 dV is 
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however, we notice that our 



121 ]. when they adopted the 



the astrophysical factor for a single subhalo with mass m. All subhalo parameters are taken 
identical to those used to account for the diffuse 7 rays. 

We then calculate, in the same propagation model, the p spectrum at Earth by incorpo- 
rating the source term of Eq. (j3j) in GALPROP. Fig. [H] shows the background, the signal 
and the total p fluxes. The background antiproton flux s lower than the data; this has 
been discussed in jgl This is another hint at the necessity of an exotic signals besides 
the ordinary CR secondaries. The p flux in our model is about one order of magnitude 
greater than the measured data at energies lower than 1 GeV. 
prediction is a few times lower than that by Bergstrom et al. 
median set of propagation parameters, which are similar to the propagation parameters in 
our conventional model. The difference may be due to two reasons: first we adopt a differ- 
ent propagation model; second — and maybe more important — , we do not adopt a universal 
"boost factor" for the diffuse 7 rays from DMA as de Boer and Bergstrom et al. did. The 
enhancement by subhalos tends to boost the 7 ray and p fluxes at large radii, as shown in 
Fig. [21 Therefore subhalos boost 7 rays more than p since only these p produced within the 
Galactic diffusion region contribute to the p flux on Earth. 

However, these effects are not enough to give a p flux consistent with the data. This is 
related to the presence of two DM rings near the solar system, that enhance the antiproton 
flux greatly as they are strong p sources. We have to resort to a new propagation model to 
suppress the p further and give a consistent description of all data. 

IV. A NEW PROPAGATION MODEL 

Inheriting the advantages of our conventional model in the last section, a new propagation 
model is intended to be built to account for the Galactic diffuse 7 rays, and at the same 
time the p flux. Degeneracies exist between the propagation parameters, mainly between 
the diffusion coefficient D(p) and the height of the diffusion region z^. Different sets of 



parameters can all explain the CR data but lead to very different signals from DMA [52 ]. 
If we adopt smaller z^ we can adjust D(p) at the same time to give similar prediction of 
the CR data. However adopting smaller Zh leads to smaller p flux from DMA since only p 
sources within the diffusion region contributes to the flux on the Earth. The degeneracy 
leads to about one order of magnitude uncertainties of the p flux when adopting different 
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sets of propagation parameters in Bergstrom et al. 12] . Therefore it is straightforward to 
consider a smaller height of the diffusion region to suppress the DMA signal. 

The diffusion halo height is determined by fitting the CR data. It is found that the 
average gas density CRs crossed during their travel to the Earth is about 0.2 atoms/cm 3 , 
which is significantly lower than the average gas density in disk of about 1 atom/cm 3 . This 
is explained by the possibility that CRs are confined in a larger diffusion region than the 
gaseous disk, spending a longer time outside the gaseous disk. The height of the halo 
represents the volume of this diffusion region. In [53J, Strong and Moskalenko derived 
Zh = 3 — 7 kpc for the four types of radioactive nuclei isotopes of 54 Mn/Mn, 36 C1/C1, 
26 A1/ 27 A1 and 10 Be/ 9 Be, while the data of 10 Be/ 9 Be favored a smaller halo of 1.5 — 6 kpc. 
Maurin et al. found that several settings of propagation parameters with halo height Zh 



ranging from 1 kpc to 15 kpc could well fit the observational B/C ratio data [54|, |55 |. 
Certainly the diffusion halo height should be at least as large as the scale height of the 
ionized component of the interstellar gas, ~ 700 pc (5^]. The halo height is also constrained 
by the diffuse 7 ray emission at the intermediate latitude, since too small halo height will 
lead to too low 7 ray flux at at these latitudes. 

Since the gas is not smoothly distributed in the disk, CRs may travel in low density 
regions in the disk and may not diffuse to so large a region as usually considered. As shown 



m 



571 ]. the gas density and the strength of magnetic fields are higher than average within 



molecular clouds: the molecular clouds can act as magnetic mirrors to reflect and confine the 
charged CR particles. In this scenario CRs may travel in a low gas density region rather than 
in a region with the average density. Nevertheless, the effect of molecular cloud reflection 
is hard to quantify. We assume the gas density that CRs crossed is lowered compared 
with the average gas density in the disk by a constant factor, which should be of the order 
of 1 since the conventional diffusion model has been very successful in describing the CR 
transportation. We find that a reduction factor ~1.5 and Zh = 1.5 kpc can reproduce all 
the data very well, which will be shown below. It should be noted that the reduction of the 
gas density crossed by CRs is the key point of this propagation model. 

Taking this smaller halo height, we then adjust the diffusion coefficient and the Alfven 
speed, so that we can reproduce the secondary to primary ratios, e.g. the B/C ratio. In the 
present model, parameters in Eq. ([3]) are taken as D = 1.3 x 10 28 cm 2 s _1 , po = 4GV and 
5 = 0.34. The Alfven speed va is taken as 19 km s _1 . The injection spectra of nuclei and 
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electrons are given in Table [B 

The CR source distribution adopted in this model is shown in Fig. HI together with 



the source distribution adopted in 381, which is the same as the pulsar distribution 142], 

n n 

and that adopted in [5|, [6J , which is obtained by fitting the EGRET data with a constant 
Xco value. The vertical bars are SNR data points from [43] . The flat source distribution 
derived by fitting the EGRET diffuse 7 ray data is not in agreement with the distribution 



of SNRs 



381 ] . As for the pulsar distribution, it is so peaked that it would aggravate the 



problem of reproducing the relatively smooth EGRET flux profile along the Galactic Plane 
. In addition, the pulsar distribution is probably not sufficiently reliable to trace the 



SNRs accurately. In fact, only a proportion of about one in ~ 100 known pulsars appears 



to be convincingly associated directly with SNRs [59|. It is very interesting that in our 
new model, the source distribution consistent with SNRs data (solid line in Fig. H]) can give 
a better description of the diffuse 7 ray data with a constant CO-to-H 2 conversion factor 
Xco ~ 1-0 x 10 20 molecules cm -2 / (K km s _1 ) than the other two distribution functions. The 
radial source distribution has the same form as the previous one with a = 1.35, j3 = 2.7. 

In Fig. [7J we show the results of B/C, 10 Be/ 9 Be, spectra of protons and electrons as 
predicted by this new propagation model. The stable secondary to primary ratio B/C and the 
unstable to stable secondaries 10 Be/ 9 Be show that our model reproduces the experimental 
date very well. In this model the interstellar proton spectrum is taken to be the locally 
observed spectrum since the energy loss of proton is negligible. The interstellar electron 
spectrum has an intensity normalization at high energies different from the local one, similar 
to the electron spectrum adopted by Strong et al. in |6j but with smaller fluctuations. We 
will show below that this model can also reproduce the diffuse 7 ray and p data well. 



A. Diffuse 7-rays and p flux 

In order to give 7 ray and p fluxes consistent with experimental data, we also changed 
the DM ring parameters slightly: the inner ring is now located at R = 3.5 kpc and the outer 
ring is moved from R = 14 kpc to 16 kpc. On one hand, the position of rings are not crucial 
parameters to fit the EGRET data: indeed, a slight change will not vary the prediction of 
diffuse 7 rays much. On the other hand, the rings are mainly helpful to explain the diffuse 7 
rays at intermediate latitudes; there are not enough to account for the GeV excess elsewhere, 
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FIG. 7: The CR spectra in our new model. Curves and data are the same as in Fig. O 



such as in regions D and E (if only the smooth DM component is included). However, our 
consideration of subhalos helps to enhance 7 ray emissions at large radii as shown in Fig. [3 
Therefore, we do not need so large 7 ray emissions from the two rings to contribute to the 
intermediate latitudes as done in the de Boer model, and we can move them slightly far 
away from the solar system. It is interesting to note that the analysis of the HI gas flaring 
by Kalberla et al. also favored a DM ring located at a large radius of ~ 17.5 kpc 60]. 
We have checked that these soft modulation of ring parameters did not change the rotation 
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FIG. 8: The rotation curve in our new propagation model. 

curve significantly. Fig. [8] shows the rotation curve for the present model. The contributions 
from the DM smooth halo, the two DM rings, and the bulge and disk are included. It can 
be seen that the rotation curve is consistent with data. 

Adding the background diffuse 7 rays from CRs and those from DMA directly, we find 
the calculated diffuse 7 rays are well consistent with observations. The diffuse 7 ray spectra 
in the six sky regions are shown in Fig. [9j The prediction in the new propagation model is 
in good agreement with EGRET data. 

Figs. [10], [UJ and [12] display the diffuse 7 longitudinal and latitudinal profiles in our 
present model. The line styles in these figures are the same as those given in Figs. [5] and 
[9] representing contributions from tt° decay, IC, bremsstrahlung and EGRB respectively. 
The solid red line is the total contribution from the background 7 rays (solid blue) and 
DMA signals. From these figures, it is obvious that the DMA component is essential to 
explain the profiles above 500 ~ 1000 MeV. We notice that the longitudinal profiles at low 
latitude (|6| < 5°) and the latitudinal profiles in the outer Galaxy with 30° < I < 330° are 
in fairly good agreement with the EGRET data. However, for the latitudinal profiles in the 
inner Galaxy (330° < / < 30°), we also find that the intensities at intermediate latitudes 
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FIG. 9: Diffuse gamma ray spectra in the six sky regions predicted by the new propagation model 
(top row, regions A, B, middle C, D, bottom E, F). The various contributions are ir° decay, inverse 
Compton, bremsstrahlung, EGRB and DMA (dark red curve). 
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FIG. 10: Longitude profiles at low latitudes (|6| < 5.5) in the new propagation model, compared 
with EGRET data in 10 energy ranges from 30 MeV to lOGeV. 
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FIG. 11: Latitude profiles in the inner Galaxy |/| < 30° in the new propagation model, compared 
with EGRET data in 10 energy ranges from 30 MeV to lOGeV. 
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FIG. 12: Latitude profiles in the longitude ranges 30° < I < 330° in the new propagation model, 
compared with EGRET data in 10 energy ranges from 30 MeV to lOGeV. 
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FIG. 13: Flux of antiprotons after solar modulation in the new propagation model. Lines and data 
are the same as in Fig. [6j 

10° ^ b 30° are lower than measurements at energies from 500 MeV to 4000 MeV. A 
similar excess is present in the "optimized model" {(J, in which it is pointed out that this 
may be related to an underestimate of the ISRF in the Galactic halo and that a factor of 
~ 2 uncertainty on the ISRF is quite possible due to the complexity of its calculation 
Furthermore, the smaller Zh adopted here shallows the distribution of electrons, which may 
also contribute to this discrepancy. 

Finally, the antiproton flux is given in Fig. [13l Below ~ 8 GeV, the p flux from DMA 
dominates the CR secondary p. The total p flux is a bit higher than the best fit value of 
the experimental data 2 . However due to large errors of the present experimental data, the 
prediction is still consistent with data within la. Forthcoming high precision p measurements 



2 Here we take the solar modulation potential <!> = 900 MV. In fact, the modulation parameter is not a free 
parameter for different CR species, but depends on the solar activities. The measurements of the p flux 
were taken from the solar maximum (BESS00, BESS02) to the solar minimum (BESS95). 
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by PAMELA 61[ and AMS02 62] will be helpful to determine if the present model is viable. 
It should be noted that our model has the potential to further suppress the p flux, by 
changing the rings position and the diffusion height z^. 



In 



631 ]. de Boer et al. introduced an anisotropic propagation model to greatly suppress 



the p flux. In their model, only the p from DMA within the gaseous disk is confined by 
magnetic field and contributes to the p flux at Earth. The p from DMA above the gaseous 
disk is blown away and has no contribution to the local p flux. In order to reproduce the 
B/C data, they introduced a grammage parameter c = 12 so that the secondaries and the 
resident time are increased by this factor, as a result of the molecular clouds confinement 
of charged particles. However, we think it may be hard to reproduce the diffuse 7 ray 
data at intermediate and high latitudes in this model, since the CRs above the gas disk 
are quickly blown away and produce very low 7 ray emissivity. Our model gives consistent 
Galactic diffuse 7 rays and p flux with experimental data simultaneously without drastic 
modifications of the GALPROP model. 



V. SUMMARY AND DISCUSSION 

In this work we propose to solve the "GeV excess" problem of the Galactic diffuse 7 rays 
by developing a new propagation model which includes contributions from DMA. We have 
shown that this propagation model can well reproduce the B/C, 10 Be/ 9 Be data and spectra 
of protons and electrons. The Galactic diffuse 7 ray spectrum at different sky regions and 
its profile as a function of longitude and latitude are also shown to describe the EGRET 
data very well, if the DMA contribution is included. The p flux in this model is consistent 
with experimental data within la. Compared with previous works jf], Q], our model does 
not introduce a normalization of the interstellar proton intensity different from the local one 
(as it should be universal due to the negligible energy loss of protons). Furthermore, neither 
the "boost factors" to the DMA signals nor the arbitrary renormalization of the Galactic 
7-ray background are needed in our model to explain the EGRET data. 

In our model, the p flux from DMA is suppressed by the following changes compared 
with the de Boer model [8|: 1) the smaller Zh helps to suppress the p flux from the smooth 
component of DM. The average gas density CRs crossed may be smaller than the average 
gas density in the disk, as shown in {5^ ■ This fact tends to favor a smaller height of the CR 
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diffusion region. 2) We do not adopt a universal "boost factor" for 7 rays and p. Since the 
enhancement by subhalos is larger at large radii, the boost of p at the solar neighborhood is 
smaller than the boost of 7 rays. 3) The ring parameters are slightly adjusted, which does 
not change the 7 ray profile while greatly suppressing the p flux. This is because the distance 
dependence of the p propagation is steeper (exponential decrease) than r~ 2 of 7-rays [64]. 

A potential problem of this model is related to the p flux, which is consistent with data, 
but cannot best fit the data. Adjusting the propagation or DM parameters (such as the 
rings) can further lower the p flux. However, more fine-tuning is required to fit the 7-ray 
and rotation curve data. Note a new explanation of the "GeV excess" was given as an 
instrumental biases contaminating the EGRET data 65] . However, whether this conjecture 
is correct or not cannot be confirmed at the moment. Forthcoming precise observations, 



such as GLAST |66j, PAMELA 



61| and AMS02 62], will be decisive to validate or disprove 



this model. At present, because of the fundamental importance of the DM problem, we 
think that any possible implication of DM signals deserves a serious treatment. 
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